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Abstract 

Searches for stable, hadronizing scalar quarks and gluinos are performed using 
the data collected with the ALEPH detector at LEP. Gluon splitting into a gluino 
or a squark pair is searched for at centre-of-mass energies around the Z resonance, in 
the e + e~ — > qqgg and qqqq processes. Stable squark pair production, and stop pair 
production with subsequent decays into a stable gluino, t — ► eg, are also directly 
searched for at centre-of-mass energies from 183 to 209 GeV. 

Altogether, stable hadronizing stop (sbottom) quarks are excluded up to masses 
of 95 (92) GeV/c 2 , and stable hadronizing gluinos are excluded up to 26.9 GeV/ c 2 , 
at 95% confidence level. In the framework of R-parity-conserving supersymmetric 
models in which the gluino and the stop quark are the two lightest supersymmetric 
particles, a 95% C.L. lower limit of 80 GeV/c 2 is set on the stop quark mass. 
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1 Introduction 



A search for squark production in e + e~ collisions is relevant at LEP energies because 
the stop quark and, to a lesser extent, the sbottom quark, could well be the lightest 
super symmetric partner of all standard model fermions [1]. 

Searches for squarks have already been performed by ALEPH [2, 3] in the framework 
of the MSSM, the minimal supersymmetric extension of the standard model [4], with R- 
parity conservation and the assumption that the lightest supersymmetric particle (LSP) 
is the lightest neutralino or the sneutrino. These searches yielded an absolute lower limit 
on the stop mass of 63 GeV/c 2 at the 95% confidence level. 

This absolute lower limit does not apply if the LSP is either the gluino or the 
squark itself. Supersymmetric models in which the gluino is the LSP are reviewed in 
Ref. [5]. Squarks as LSP's are cosmologically disfavoured because of their nonzero electric 
charge [6], but could be sufficiently stable to behave like the LSP in the LEP detectors. 
Scenarios in which the gluino is the next-to-lightest supersymmetric particle (NLSP) and 
decays into a b quark and a long-lived sbottom (with mass 2 to 5GeV/c 2 ) have been 
proposed to explain excesses in bb production at hadron colliders [7]. 

In this paper, it is assumed that the LSP is either the gluino or a squark, in the context 
of the MSSM with R-parity conservation. Under this assumption, the NLSP would decay 
with a 100% branching fraction into the LSP, e.g., t — > eg, b — > bg or g — > bb. (Stop 
decays into tg and gluino decays into tt can only happen for stop and gluino masses larger 
than the top quark mass, which is beyond the reach of LEP.) This coloured LSP is stable 
and hadronizes into stable, colourless, charged and neutral bound states. These bound 
states (e.g., gg, gqq, gqqq, qq, qqq) are called R-hadrons, in which "R" refers to the fact 
that they carry one unit of R-parity [8]. 

An important consequence is that the missing energy signature might not 
automatically be present, in contrast to other searches for supersymmetric particles 
in which the LSP is a weakly- interacting particle, i.e., a neutralino or a sneutrino. 
However, when the LSP is sufficiently heavy, the maximum energy available for R-hadronic 
interactions turns out to be quite small. Indeed, the centre-of-mass energy E scat of the 
R-hadronic interaction with ordinary matter, made of nucleons of mass wn, is very close 
to the mass mg of the R-hadron, almost independently of the R-hadron energy 
i?scat = + + 2777n-Er. It leaves little energy (E scaXi — — m N ) for each R-hadronic 
interaction in the calorimeters. For example, for m R = 50 GeV/ c 2 and E R = 90 GeV, this 
energy amounts to 800 MeV, to be compared to ~ 12 GeV for a pion of the same energy. 

Because of this reduced hadronic interaction, a large missing energy signal is expected 
from heavy neutral R-hadrons. Heavy charged R-hadrons are expected to interact mostly 
electromagnetically like heavy muons. Searches for missing energy and for stable heavy 
charged particles are therefore well suited to look for stable squarks and gluinos. 

In view of covering all possible configurations in the plane (rag, raq) with a squark or a 
gluino LSP (or, equivalently, with a long-lived squark or gluino), the processes investigated 
in this paper are 

1. the e + e~ — > qqgg process, with a gluon splitting into a pair of stable gluinos; 
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2. the e + e — ■> qqqq process, with a gluon splitting into a pair of stable squarks; 

3. the pair production of stable squarks, e + e~ — > tt and bb; 

4. the stop pair production with decays into stable gluinos, e + e~ — > tt — > cgcg. 

The first two processes were searched for using the data collected by ALEPH at LEP 1, 
at centre-of-mass energies around the Z resonance. These data correspond to about 4.5 
million hadronic Z decays. The data collected at LEP 2 were used to analyse the last 
two processes. The integrated luminosities and centre-of-mass energies of these data are 
indicated in Table I. 

Table 1: Integrated luminosities, centre-of-mass energy ranges and mean centre-of-mass energy values 
for the data collected during the years 1997-2000. 



Year 


Luminosity [pb 1 




Energy range [GeV] 


(V5> [GeV] 


2000 


9.4 


207-209 


208.0 




122.6 


206-207 


206.6 




75.3 


204-206 


205.2 


1999 


42.0 




201.6 




86.2 




199.5 




79.8 




195.5 




28.9 




191.6 


1998 


173.6 




188.6 


1997 


56.8 




182.7 



This paper is organized as follows. The detector is briefly described in Section 2. 
The simulation of the signal final states is discussed in detail in Section 3, and a brief 
account of the simulation of the standard model backgrounds is given. The e + e~ — > qqgg 
and qqqq searches at LEP 1 are presented in Sections 4 and 5, followed by the search 
for stable squark pair production at LEP 2 in Section 6, and by the search for decaying 
stop quarks in Section 7. The result of the combination of all these analyses is given in 
Section 8. 

2 The ALEPH detector 

A detailed description of the ALEPH detector and of its performance can be found in 
Refs. [9, 10]. Only a summary is given here. 

Charged particles are detected in the central part, consisting of a precision silicon 
vertex detector, a cylindrical drift chamber and a large time projection chamber, which 
measure altogether up to 31 space points along the charged particle trajectories. The time 
projection chamber also provides 359 measurements (338 from wires and 21 from pads) of 
the specific energy loss by ionization dE/dx. A 1.5 T axial magnetic field is provided by a 
superconducting solenoidal coil. Charged particle transverse momenta are reconstructed 
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with a 1/pt resolution of 6 x 1CT 4 05 x 10~ 3 /p T {p T in GeV/c). In the following, good 
tracks are defined as charged particle tracks reconstructed with at least four hits in the 
time projection chamber, originating from within a cylinder of length 20 cm and radius 

2 cm coaxial with the beam and centred at the nominal collision point, and with a polar 
angle with respect to the beam such that | cos#| < 0.95. 

Electrons and photons are identified by the characteristic longitudinal and transverse 
developments of the associated showers in the electromagnetic calorimeter, a 22 radiation 
length thick sandwich of lead planes and proportional wire chambers with fine read-out 
segmentation. The relative energy resolution achieved is 0.18/ y/E (E in GeV) for isolated 
electrons and photons. 

Muons are identified by their characteristic penetration pattern in the hadron 
calorimeter, a 1.5 m thick yoke interleaved with 23 layers of streamer tubes, together with 
two surrounding double-layers of muon chambers. In association with the electromagnetic 
calorimeter, the hadron calorimeter also provides a measurement of the hadronic energy 
with a relative resolution of 0.85 /\/E (E in GeV). 

The total visible energy, and therefore the missing energy, is measured with an energy- 
flow reconstruction algorithm which combines all the above measurements, supplemented 
by the energy detected down to 34mrad from the beam axis (24mrad at LEP 1) by 
two additional electromagnetic calorimeters, used for the luminosity determination. The 
relative resolution on the total visible energy is 0.&0/\fE (E in GeV) for high-multiplicity 
final states. This algorithm also provides a list of reconstructed objects, classified 
as charged particles, photons and neutral hadrons, called energy-flow particles in the 
following, and used to determine the event characteristics for the selections presented in 
this paper. 

3 Monte Carlo Simulation 

3.1 Signal simulation 
3.1.1 Production processes 

The simulation of squark pair production at LEP 2 energies was performed with the 
PYTHIA event generator [11]. For the simulation of gluon splitting into a gluino pair at 
the Z resonance, the program mentioned in Ref. [5] was used, modified to include initial- 
state radiation as described in Ref. [12] and to generate events as in Ref. [13]. The QCD 
leading-order production cross section for e + e~ qqgg was determined numerically with 
this program, and cross-checked with the analytical predictions of Refs. [14, 15], after 
proper quark to gluino colour-factor modification. These analytical formulae also allowed 
the e + e~ — > qqqq cross section to be determined, by appropriately modifying the quark to 
a squark phase space factor. The calculation of the resummed QCD next-to-leading-order 
cross section for gluon splitting into a gluino or a squark pair was performed following 
the prescriptions of Ref. [15], with the running of the strong coupling constant corrected 
for the presence of a light gluino or a light squark [16]. The resulting Z partial widths are 
displayed in Fig. 1 as a function of the gluino and the squark mass. 
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Figure 1: The Z — > qqgg and Z — > qqqq partial widths as a function of the squark and the gluino 
mass, computed at the leading order (dashed curves) and resummed at the next-to-leading order (full 
curves). Also shown are the corresponding 95% C.L. lower limits on the squark and gluino masses when 
the contribution to the Z width exceeds 3.9 MeV (Ref. [17] and Section 4). 



3.1.2 R-hadron formation 

The hadronization of the (s)partons generated as described above was done with PYTHIA. 
Because squark and gluino hadronization is not available by default, PYTHIA had to be 
extended with a few dedicated routines [18]. The following issues are addressed therein. 

Stable gluino and squark hadronization Stable gluinos and squarks are allowed to 
radiate gluons [19] as any other coloured particles. The fragmentation is handled with a 
Peterson function [20] with a parameter extrapolated [21] from its value for b quarks: 

£M = _^_. (1) 

However, a gluino is attached to two, rather than one, string pieces. In this case, the 
fragmentation is applied to each of the two pieces, thereby giving a larger energy loss 
than for a (s) quark. 

Both R-meson and R-baryon production are possible. All R-hadrons made of a given 
squark or gluino are assumed to be equal in mass, with an electric charge of —1, or 1, 
determined by the quark and squark content. The relative fractions of charged and neutral 
R-hadrons are then obtained by simple statistics, and are typically half and half, except 
for gluino R-hadrons, for which the additional possibility to form a gluino-gluon bound 
state enhances the fraction of neutral R-hadrons by an unknown amount. 
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In the processes studied in this paper, R-hadrons are produced in pairs, so that double 
neutral (R°R°), mixed (R°R ± ) or double charged (R ± R ± , R ± R =F ) final states are possible. 
There is no significant correlation between the electric charges of the two gluino R-hadrons. 
Charge conservation is ensured by the presence of additional fragmentation particles in 
the final state. 

Stop-hadron decay The stop decay to eg proceeds through loop diagrams or through 
unlikely FCNC tree-level couplings [22], which leads to a lifetime larger than the typical 
hadronization time. The stop quark therefore hadronizes into a stop R-hadron before 
decaying. The decay is described in the framework of the spectator model [23] , in which 
the bound state quarks act as spectators to the decaying stop. 

3.1.3 R-hadronic interaction in the calorimeters 

A major issue of the simulation is the treatment of the interaction of R-hadrons in the 
detector. Extensive studies have been done on this subject [5, 24]. The simple approach 
chosen here to simulate R-hadronic interactions, in analogy with Ref. [3], is to treat R- 
hadrons as heavy pions in the GHEISHA package (in which the hadron mass is properly 
accounted for), with the modification that the low energy pion-nucleon resonances are 
removed. 

3.1.4 Free parameters 

Several unknown parameters, listed below, are needed to fully specify the R-hadron 
phenomenology, and may be the source of sizeable systematic uncertainties. 

• The probability P gg to form a gluon-gluino bound state. This quantity is essentially 
unknown and can vary between and 1. For P gg = 0, the fractions of double 
neutral, mixed and double charged final states are approximately 25%, 50%, and 
25%, respectively. In the configuration in which P gg = 1 (unlikely for gluino masses 
in excess of 25GeV/c 2 [5, 24]), only purely neutral final states are produced. 

• The effective spectator mass M^ ec of squark and gluino R-hadrons. It 
roughly corresponds to the difference between the constituent mass and the 
current algebra mass. The value obtained from hadron mass spectroscopy is 
<L = 0.5i^GeV/c 2 [25]. 

• The gluon constituent mass M|. Its value is estimated from glueball searches [26] 
to be M| = 0.7±o;4 GeV/c 2 , to be compared to the quark constituent mass, typically 
twice smaller. 

• The total cross section <7rn of R-hadron-nucleon scattering. It is assumed to be equal 
to the cross section of pion-nucleon scattering because (i) at high momentum, 
the interaction cross section of any hadron on a nucleon is proportional to the 
sum of the individual valence parton cross sections, and (ii) heavy partons (with 
a mass above a couple of GeV/c 2 ) do no interact significantly. The sum therefore 
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runs only over the standard partons, in which each quark accounts for one unit of 
cross section, and each colour octet gluon for 9/4 units. For gluinos, the R- meson 
(gqq) hadronic cross section is thus equal to that of a pion, the gg bound state 
cross section is 9/4/(1 + 1) = 9/8 larger, and the R-baryon (gqqq) cross section is 
(1 + 1 + 1)/(1 + 1) = 3/2 larger. For squark-R-mesons (qq), the same argument 
yields a ratio of 1/2 while, for squark-R-baryons (qqq), the ratio is equal to unity. 
Altogether, it is therefore reasonable to use an average R-hadronic interaction cross 
section identical to that of pions, with an uncertainty of ±50%. 

3.1.5 Simulated signal samples 

In order to design the signal selection criteria, hundreds of simulated signal samples of 1000 
to 4000 events each were generated for each of the four processes mentioned in Section 1. 
Gluon splitting to gluino or squark pairs was simulated for gluino/squark masses between 
and 40GeV/c 2 , with a typical step of 2GeV/c 2 . For squark pair production, samples 
with squark and gluino masses ranging from to 90GeV/c 2 with a typical 5GeV/c 2 step 
were generated, with two values of the squark mixing angle: (i) 0° and (ii) the value for 
which the squark coupling to the Z vanishes, i.e., 56° for stops and 68° for sbottoms. 

3.2 Background simulation 

Simulated samples for all relevant standard model background processes were generated 
both at LEP 1 and LEP 2 energies. Bhabha scattering was simulated with the BHWIDE 
generator [27]. The K0RALZ package [28] was used for the other difermion processes. 
Two-photon interactions were simulated with PH0T02 [29]. The W + W~ production was 
simulated with K0RALW [30], the production of Wez/ with GRC4F [31] and the Zee, ZZ and 
Zi/P final states with PYTHIA [32]. 

4 Search for e + e _ — ► qqgg with LEP 1 data 

At LEP1, events resulting from gluon radiation off a qq pair, with subsequent gluon 
splitting into two stable gluinos, are expected to show as a pair of acoplanar jets, 
accompanied by two stable R-hadrons. 

When the gluino mass is small, the e + e~ — > qqgg production cross section is large 
enough to sizeably contribute to the Z hadronic width. The accurate electroweak 
measurements at LEP and SLC allow a model-independent upper limit of 3.9 MeV to 
be set on the Z width for purely hadronic final states [17]. All gluino masses below 
6.3GeV/c 2 are therefore excluded at the 95% confidence level [17], as can be seen in 
Fig. 1, irrespective of the gluino decay and hadronization mechanisms. For larger masses, 
the final state topology, and therefore the search strategy, depends on the electric charges 
of the R-hadrons. 
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4.1 Search for two neutral R-hadrons 



When the two R-hadrons are neutral, the acoplanar jet pair is accompanied by missing 
energy arising from the specific interaction of massive neutral R-hadrons with the detector. 
The acoplanar jet selection developed for the Hi/z/ search [33] may therefore be used. Its 
efficiency reaches ~ 9% for mg = 25 GeV/c 2 , as is visible from the dashed curve in Fig. 2. 
This analysis was optimized for the case of a heavy Higgs boson, leading to non- violently 
acoplanar hadronic jets and a moderate, but isolated, missing energy, which makes it 
inadequate for gluino masses above 30 and below 15GeV/c 2 . 
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Figure 2: The qqgg selection efficiency as a function of the gluino mass, for the Hvv selection (dashed 
curve), the additional efficiency of the large-mass selection (dot-dashed curve), the additional efficiency 
of the small-mass selection (dotted curve), and the sum of the three (full curve). The markers with error 
bars indicate the efficiencies obtained with each of the simulated samples, and the curves are obtained 
with a polynomial fit through these points. 



4.1.1 Large gluino masses 

For gluino masses above 30GeV/c 2 , the missing energy and the jet acoplanar ity become 
so large that simpler and, in this configuration, more efficient selection criteria can be 
designed with the same topological variables as in the Yivv search. Such selection criteria 
were optimized for a search for the stop quark [34], and are summarized below. 

Only events with at least four good tracks are considered. Most of the qq background 
is rejected by the requirement that the visible mass be smaller than 50%-y/i. Events with 
a small angle 9 T of the thrust axis with respect to the beam (cos^t > 0.9) are eliminated, 
and so are events with energy detected within 12° from the beam axis, to avoid events 
in which energy escapes down the beam pipe. The bulk of two-photon interactions is 
rejected by requiring that the total momentum transverse to the beam exceed 5%y/s. 
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For events with at least eight good tracks (mostly qq events) , the acollinearity and the 
acoplanarity angles are required to be smaller than 135° and 150°, respectively. Three-jet 
events with one energetic neutrino from heavy-quark semi-leptonic decays are rejected by 
requiring that the total momentum transverse to the beam remain smaller than 30 GeV/ c. 
To reject three-jet events with two such neutrinos, the aplanarity (i.e., the sum of the 
three jet-jet angles when the event is forced to form three jets) has to be less than 350°. 

In events with less than eight good tracks, the thrust axis angle with respect to the 
beam axis has to be larger than 45°. For monojet events, i.e., events in which one 
hemisphere (with respect to a plane perpendicular to the thrust axis) is empty, it is 
further required that the visible mass be in excess of 5 GeV/ c 2 and the thrust smaller than 
0.95. These cuts eliminate the two-photon interactions and the e + e~ — > r + r~ background. 
For non-monojet events, each hemisphere must contain at least two good tracks and the 
acoplanarity angle has to be less than 150°, which rejects the remaining r + r" background. 

The additional efficiency brought by this large-mass selection exceeds 10% for mg 
above 30 GeV/ c 2 , as can be seen from the dot-dashed curve in Fig. 2. 

4.1.2 Small gluino masses 

For gluino masses below 15 GeV/ c 2 , the R-hadrons tend to interact more in the 
calorimeters, and a third jet, mostly formed with neutral hadronic energy, develops around 
the still large missing energy. Requirements like missing energy isolation or substantial 
aplanarity therefore become inefficient, and must be replaced by more selective criteria 
aimed at this specific three-jet topology. 

Only hadronic events (more than five good tracks carrying more than 10% of the 
centre-of-mass energy) with a visible mass smaller than 65 GeV/ c 2 are considered here. 
As in the Hz/z/ selection, the acollinearity and the acoplanarity angles are required to be 
less than 165° and 175°, respectively. Similarly, events with less than 75% of the visible 
energy above 30° from the beam axis are eliminated, and so are the events with any 
activity below 12°. Two-photon interactions are rejected by requiring that the missing 
transverse momentum exceed h%\fs for events with a visible mass below 20 GeV/ c 2 , and 
that the missing momentum point more than 35° away from the beam axis. The latter 
cut also rejects qq events with an undetected, energetic, initial-state-radiation photon. 

The remaining events are then forced to form three jets with the Durham algorithm. 
To ensure a three-jet topology, the values where the transitions from three to two jets 
(1/23) and to four jets (1/34) occur are required to exceed 0.01 and to be less than 0.04, 
respectively. 

To select signal-like events, the "R-hadron jet" candidate, i.e., the jet with the smallest 
amount of charged energy, is required to contain mostly neutral hadrons, by imposing 
that its charged energy be smaller than 20% of its neutral hadronic energy. (The latter 
is computed as the energy sum of all neutral energy-flow particles in the jet with a 
hadron calorimeter component.) Because the missing energy is also expected to point in 
the direction of the "R-hadron jet" candidate, it is finally required that the amount of 
charged energy in a cone of 35° around the missing momentum direction be smaller than 
3 GeV. As can be inferred from Fig. 3, these last two criteria cut the expected background 
by a factor larger than a thousand, but only reduce the signal efficiency by a modest 30%. 
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Figure 3: Distributions of (a) the ratio of the "R-hadron jet" charged energy to its neutral hadronic 
energy; and (b) the charged energy in a 35° cone around the missing momentum, before the corresponding 
cuts are applied, for the data (triangles with error bars), the background simulation (light-shaded 
histogram) and the signal with m g — 10 GeV/ c 2 (dark-shaded histogram, arbitrary normalization). 

The additional efficiency brought by this last selection is shown by the dotted curve 
in Fig. 2. As anticipated, the selection is well suited for masses around 15 GeV/ c 2 and 
below. Although not specifically designed for large masses, the selection also contributes 
substantially all the way up to 40 GeV/ c 2 . 

When the three searches are combined, the qqR°R° selection efficiency exceeds 15% 
for masses between 10 and 35 GeV/c 2 , and culminates at about 35% for mg = 27 GeV/ c 2 . 
It decreases fast for small gluino masses but it does not completely vanish, even for 
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massless gluinos. Indeed, when two energetic light neutral R-hadrons are produced, they 
are expected to be detected as normal neutral hadrons in the hadron calorimeter. Due 
to calorimeter resolution effects, this large neutral hadronic energy occasionally causes 
sizeable missing energy and acoplanarity. 

No events were selected in the data in either of the three selections, with an expected 
background of less than two events at 95% confidence level. 

4.2 Search for charged R-hadrons 

If the probability Pg g to form a gluon-gluino bound state is not 100%, a relevant fraction 
of the final states contains one or two charged R-hadrons (up to 75% for Pg g = 0). These 
heavy stable charged particles can be identified with their isolation from the rest of the 
event, as well as their unusual energy loss by ionization in the time projection chamber. 

Events with at least five good tracks are considered in this search. Heavy stable 
charged particle track candidates must be reconstructed with at least half of the 338 
wires and the 21 pads of the time projection chamber, must have a momentum larger 
than 2.5GeV/c, and must be isolated from the other good tracks of the event by more 
than 18.2°. To ensure a reliable momentum determination, the track-fit-x 2 probability is 
required to exceed 1%. The energy loss by ionization is then required to be at least three 
standard deviations away from that expected for light stable charged particles (electrons, 
muons, pions, kaons or protons) of the same momentum. After this preselection, a total of 
364 events with at least one heavy stable charged particle track candidate is selected from 
the data, in agreement with the 363 ± 16 events expected. The distribution of N a , the 
number of standard deviations with respect to the ionization expected in the most-likely 
light-particle hypothesis, is displayed in Fig. 4a, and that of the cosine of the isolation 
angle in Fig. 4b. A reasonable agreement is observed between the data and the simulation 
in both variables. In particular, the dE/dx distribution of very ionizing particles (mostly 
alpha particles) is well reproduced by the simulation. 

The distributions of N a as a function of the track candidate momentum are shown in 
Fig. 4c for the data, and in Fig. 4d for signal events with different gluino masses. Events 
with two heavy stable charged particle track candidates are kept as qqR ± R ± candidates. 
Events with only one such track, i.e., qqR ± R° candidates, are required to have a visible 
mass smaller than 75GeV/c 2 , to account for the presence of a neutral R-hadron. The 
track candidate must also satisfy one of the two following tightened identification criteria 
in order to reject the remaining background: 

(i) either the track is isolated from the other good tracks of the event by more than 
25.8° and its fit-x 2 probability is in excess of 10%, in which case its energy loss 
by ionization is required to be at least four standard deviations away from that 
expected for light stable charged particles of the same momentum; 

(ii) or its momentum is larger than 5 GeV/c, in which case the energy loss by ionization 
is required to be smaller than expected for light charged particles. 
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Figure 4: Distributions of (a) \N a \, the number of standard deviations with respect to the expected 
dE/dx for the most-likely light-particle hypothesis; and (b) the cosine of the isolation angle, for the 
heavy-stable-charged-particle track candidates after the preselection, in the data (triangles with error 
bars) and the background simulation (histogram). Distributions of N a as a function of the charged 
particle momentum (c) in the data and (d) in the signal simulation with mg = 4, 10 and 25GeV/c 2 . 



4.3 Combined selection efficiency 

The efficiencies of the combination of all (charged and neutral) R-hadron selections, 
weighted by the expected production fractions, are displayed in Fig. 5 as a function 
of the gluino mass for P gg = 0%, in each of the three different final states, qqR°R°, 
qqR ± R° and qqR ± R ± . The total efficiency is maximal for gluino masses in the vicinity of 
27GeV/c 2 , where it exceeds 50% (to be compared to 35% when P gg = 100%). No events 
were selected in the data by the charged hadron selection, with 0.7 events expected from 
standard model backgrounds, mostly from e + e~ — > qq. 
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Figure 5: The qqgg selection efficiency as a function of the giuino mass, when Pg g = 0%, for the 
qqR°R° final state (dashed curve), the qqR ± R° final state (dot-dashed curve), the qqR ± R ± final state 
(dotted curve), and the sum of the three (full curve). The markers with error bars indicate the efficiencies 
obtained with each of the simulated samples, and the curves are obtained with a polynomial fit through 
these points. 



4.4 Systematic uncertainties 

Potentially large systematic uncertainties may arise from several sources and are listed 
below. Their relative effects on the number of signal events expected are given here for a 
giuino mass of 27GeV/c 2 , but are similar for all masses. 

• The hadronization is simulated with parton shower evolution, but the parton shower 
parameters are mostly unknown in the presence of gluinos. A variation of these 
parameters by ±100% yields a ±3% variation of the selection efficiency 

• The R-hadronic interaction cross section uncertainty of ±50% changes the selection 
efficiency by ±9% for Pg g = 1 and by ±3% for Pg g = 0. 

• The gluon constituent mass uncertainty yields an efficiency variation of ±2% for 
Pg g = 1. The other giuino hadronization parameters of Section 3 have no sizeable 
effects. 

• The uncertainty on as(mz) = 0.1183 ± 0.0020 turns into a variation of ±4% of the 
production cross section. 

• The QCD next-to-leading-order corrections to the cross section are at the percent 
level for giuino masses around 25GeV/c 2 , and are included in the present estimate. 
The next-order corrections are not expected to have any visible effect. 

The above uncertainties were taken into account according to Ref. [35] . 
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4.5 Gluino mass limit 



The number of qqgg events expected in the LEP 1 data is displayed in Fig. 6 as a function 
of the gluino mass, for Pg g varying between and 1. Altogether, when combined to the 
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Figure 6: The number of qqgg events expected to be selected in the LEP 1 data by the neutral and 
charged R-hadron searches, as a function of the gluino mass, for P gg varying between and 1 (shaded 
area). The markers with error bars (triangles for P gg — 1, stars for P gg = 0) show the individual values 
obtained from each of the simulated signal samples. Also shown is the region excluded by the Z lineshape 
measurement (hatched region), and the lower limit on the gluino-LSP mass (dot-dashed lines). 

analysis of Ref. [17], this search results in a 95% C.L. absolute lower limit of 26.9 GeV/c 2 
on the mass of a gluino LSP (Fig. 6), which substantially improves the 2 to 18GeV/c 2 
exclusion of Ref. [36] and exceeds the expectation from the study of Ref. [5] . 

5 Search for qqqq with LEP 1 data 

The same selections were applied with no modification to stable squark production 
through the e + e~ — > qqqq process. As shown in Fig. 1, the precise electroweak 
measurements allow all squark masses below 1.3 GeV/c 2 to be excluded with the method 
of Ref. [17], irrespective of the squark decay and hadronization. 

The number of events expected to be selected by the charged and neutral R-hadron 
searches described above is displayed in Fig. 7 for larger squark masses. This search 
results in a 95% C.L. lower limit of 15.7GeV/c 2 on the mass of a stable squark. The 
result was translated into upper limits on lifetimes, as presented in Fig. 8 for squarks 
and gluinos. Conservatively, the selection efficiency was assumed to vanish for neutral- 
R-hadron decays within the calorimeter volume, and for charged-R-hadron decays within 
the tracking volume. The scenario of Ref. [7] is excluded for sbottom lifetimes in excess 
of 1 ns. 
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Figure 7: The number of qqqq events expected to be selected in the LEP 1 data as a function of the 
squark mass (full curve). The triangles with error bars show the individual values obtained from each 
of the simulated signal samples. Also shown is the region excluded by the Z lineshape measurement 
(hatched region), and the lower limit on the squark-LSP mass (dot-dashed lines). 




5 10 15 20 25 30 



M squark , M gluino (GeV/c 2 ) 

Figure 8: The upper limit on the squark and gluino lifetimes as a function of their mass (full curve). 
The hatched areas are excluded at more than 95% C.L. 
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6 Search for stable squarks in e + e — ► qq at LEP 2 

If the lightest super symmetric particle is a squark or, more generally, if the lightest 
squark does not decay within the detector volume, its pair production gives rise to a 
final state with two stable squark- R-hadrons, of which at least one is charged in at least 
60% of the cases [3]. These charged R-hadrons can then be identified using the kinematic 
characteristics of squark-pair production and the large specific ionization dE / dx measured 
with the time projection chamber, as described in the previous section. 

The stable stop search of Ref. [3] was extended to lower stop masses so as to provide 
an overlap with the search of Section 5, and to sbottom pair production. The production 
cross sections for stops and sbottoms, and the corresponding upper limits derived from 
the stable stop search, are shown in Fig. 9 for a vanishing coupling to the Z. 
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Figure 9: The production cross sections for e + e~ — > tt and bb at ^fs — 209 GeV, as a function of the 
squark mass. Also shown is the observed 95%C.L. upper limit on the cross section, derived from the 
stable stop search of Ref. [3] . 

This stable heavy charged particle search therefore allows stable stop masses between 
5 and 95 GeV/ c 2 and stable sbottom masses between 7 and 92 GeV/ c 2 to be excluded 
at 95% C.L. When combined with the result of the search for e + e~ — > qqqq with LEP 1 
data presented in the previous section, stable hadronizing squarks are excluded below 
95 GeV/ c 2 (up-squarks) and 92 GeV/ c 2 (down-squarks). 

7 Search for e + e _ — > tt — > cgcg at LEP 2 

If the gluino is the LSP and the lighter stop quark the NLSP, the stop then decays into 
a gluino and a c or a u quark. The corresponding decay widths [22] are such that, for all 
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practical purposes, the stop quark is essentially stable if the mass difference AM between 
the stop-R-hadron and the gluino-R-hadron masses is smaller than the D mass, while it 
decays promptly if AM > mo, i.e., when the t — > eg decay channel is open. In the former 
case, the stable stop search of Section 6 applies here too with no modification. In the 
latter case, the final state topology is a pair of acoplanar jets from the cc hadronization 
accompanied by missing energy carried away by the mass of the two gluinos. The selection 
of this final state and its results are discussed in Section 7.1. The AM ~ m D case, in 
which the stop lifetime can take any value, is addressed in Section 7.2. 

7.1 The case of prompt stop decays 
7.1.1 Event selection 

The discriminant variables listed below were designed for generic acoplanar jet topologies: 

• the visible mass M Y - 1S and energy E vis , computed with all energy-flow particles; 

• the energy E± 2 detected within 12° from the beam axis; 

• the total energy carried by neutral hadrons, E^n, 

• the energy of the most energetic lepton, Eg, 

• the energy in a 30° half-angle cone around the most energetic lepton, Ef°; 

• the energy computed without the identified leptons, E had ; 

• the energy measured in a 30° azimuthal wedge around the direction of the missing 
transverse momentum, i? Wedge ; 

• the energy measured beyond 30° from the beam axis, E 30 ; 

• the momentum transverse to the beam axis, computed with all energy flow particles, 
Pt, without neutral hadrons, p^ NH and with good tracks only, p^ 1 ; 

• the number of good tracks, N^; 

• the acoplanarity angle $ ac0 p, between the directions of the momenta in the two 
event hemispheres, projected onto a plane perpendicular to the beam axis. Here, 
the hemispheres are defined with respect to a plane perpendicular to the thrust axis; 

• the transverse acoplanarity angle, $ a copT, defined as above, but the hemispheres are 
now defined with respect to the transverse thrust axis, i.e., the thrust axis of the 
event projected onto a plane perpendicular to the beam axis; 

• the cosine of the polar angle of the missing momentum vector, cos# m i ss ; 

• the cosine of the polar angle of the thrust axis, cos# T ; 

• the event thrust T; 
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the polar angle of the scattered electron, S cat, computed from the missing energy 
and momentum under the hypothesis that the event originates from a 77 interaction. 



A preselection was used against two-photon events, characterized by their small visible 
energy and their boost along the beam direction. To reject these events, M vis was required 
to exceed 15 GeV/c 2 , pt to be larger than 2.5%^ and to be in excess of 7. Moreover, 
the angles $ aC op and $ a co P T were both required to be smaller than 178.5°, cos^miss and 
cos#t smaller than 0.85 and £12 smaller than 0.5%\/s. 

The distribution of the visible energy after this preselection is shown in Fig. 10 for 
the data, the background expected from different sources and the signal for two different 
values of AM. For events with small AM, only a small amount of energy is available for 
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Figure 10: Distributions of the visible energy observed in the data at y/s = 200 GeV (dots with error 
bars) and expected from various background sources (shaded histograms), after the preselection cuts 
(a). The lower two plots show the expected distributions for the signal, with AM = 5GeV/c 2 (b) and 
40 GeV/ c 2 (c). For small AM, the R-hadron composition is clearly visible: double-charged-R-hadrons 
events have largest visible energy, and double-neutral-R-hadrons events have smallest visible energy. For 
large AM, the three peaks are merged because of the broader visible energy distributions. 
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the recoiling c-quark system, which results in a small particle multiplicity and small visible 
energy. In addition, the R-hadrons tend to be emitted back to back, which turns into 
large values for thrust, acoplanarity and transverse acoplanarity. These characteristics 
are very similar to those of two-photon events. In contrast, events with large AM are 
characterized by large particle multiplicity and visible energy, and the background arises 
mainly from four-fermion production and qq events. 

In fact, the distribution of most of the variables of the above list, and thus the 
relevant background sources, are strongly correlated with the value of AM. Three different 
selections were therefore developed for small, intermediate and large AM values. 

The selection criteria follow closely those used for the squark searches in the case of 
a neutralino LSP, described in Ref. [37]. However, the values of the most relevant cuts 
were re-optimized with the N$ 5 method (that minimizes the expected 95% CL upper 
limit on the signal cross section [38]) to account for the reduced missing energy, the 
larger integrated luminosity and the centre-of-mass energy increase, and to include the 
subtraction of the four-fermion and qq backgrounds. The optimized cut values for the 
three selections are given in Table 2. 

Table 2: Selection criteria for the search for R-hadrons from stop pair production in the acoplanar jet 
topology, in the three AM regions. 



Variable 


Small AM 


Intermediate AM 


Large AM 


M vis 


> 15 GeV/c 2 


> 15 GeV/c 2 


> 15 GeV/c 2 


Pt/Vs 


> 2.5% 


>4% 


> 7.5% 


N ch 


> 7 


> 11 


> 18 


E vis / y/s 


< 50% 


< 70% 


< 80% 




< 0.5% 


< 0.5% 


< 0.5% 


COS m [ ss 


> 0.8 


> 0.8 


> 0.8 


COS 9 T 


> 0.8 


> 0.8 


> 0.85 


^acop 


< 178.5° 


< 176° 


< 176° 


^acopT 


> 178.5° 


< 177° 


< 177° 


-EWedge/ 


< 12.5% 


< 12.5% 


< 12.5% 


T 


<0.96 


<0.94 


<0.92 


^scat 


> 5° 


[6°, 80°] 


[15°, 80°] 


Pt/ -Evis 




> 12.5% 


> 12.5% 




> 1.5% 








>2% 






-^had 


[10GeV,40%v^ 


< 55% v/i 


< 75% 0v /i 




< 10% 






Enh/ E vis 




<30% 




E t /y/~S 


> 20% 






Ef°/V~s 




> 1% 


> 1% 


E30/ E v i S 






> 80% 
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7.1.2 Selection efficiency 



The selection efficiencies were obtained from simulated samples (Section 3.1.5) for over 
150 different values of (mj, mg), and were then parametrized as a function of mi and mg 
with a polynomial function. For example, the efficiency of the three selections is shown 
in Fig. 11 for mi = 80GeV/c 2 as a function of the gluino mass. The AM intervals in 
which each of the three selections are to be used were again determined, as a function of 
the stop mass, with the Ng$ prescription. For example, for mi = 80 GeV/c 2 , the switching 
AM values are 16 and 33GeV/c 2 . These switching AM values tend to decrease with the 
stop mass. 
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Figure 11: The three parametrized selection efficiencies for mj = 80GeV/c 2 (full curve, large AM 
selection; dashed curve, intermediate AM selection; and dot-dashed curve, small AM selection) as a 
function of mg. The dots with error bars show the actual efficiencies obtained from the individual 
simulated samples. 



7.1.3 Systematic uncertainties 

The signal efficiencies may be affected by uncertainties in the simulation of the stop 
and gluino R-hadronization physics, uncertainties related to the detector response for R- 
hadrons and uncertainties due to the limited size of the simulated samples. Systematic 
effects from the physics assumptions were estimated by varying the free parameters of the 
simulation, as described in Section 3. 

• The gluino-gluon bound state probability was assumed to be 10%. The effects of 
a probability variation from 0% to 100% depend on the stop and gluino masses. 
For light gluinos (around 30GeV/c 2 ), the efficiency increases with Pg g , i.e., with 
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the fraction of R°R° final states, as is the case for the LEP 1 neutral R-hadron 
selections (Section 4.1). Indeed, such light R-hadrons substantially interact in the 
calorimeters, and it is only when two neutral R-hadrons are present in the final state 
that the missing energy becomes large enough to make the present selection fully 
efficient. In contrast, for larger gluino masses and small mass differences, the two 
gluinos are produced approximately back to back with similar energies, and have 
almost no interaction in the calorimeters. This configuration leads to a cancellation 
in the measured missing transverse momentum in the R°R° final state. In this case, 
the efficiency decreases when Pg g increases. A map of the corresponding uncertainty 
(up to ±20% in extreme cases) was built in the (mi, mg) plane. 

• The R-hadronic interaction cross section uncertainty of ±50% changes the selection 
efficiency by at most ±5% in the small AM region, and by a negligible amount 
otherwise. 

• The effective spectator mass uncertainty yields relative efficiency variations of ±5%, 
±3% and ±3% in the small, intermediate and large AM regions, respectively. The 
gluon constituent mass uncertainty has no visible effect. 

• The uncertainty on the Peterson fragmentation parameter eb (varied here between 
0.003 and 0.010) gives rise to an efficiency change of ±8%, irrespective of the AM 
value. 

The limited statistics of the 150 signal samples (1000 events each) and the 
parametrization of the efficiencies with a polynomial are responsible for 3% uncertainty. 
Beam- related backgrounds, which affect the determination of E i2 , were not simulated. 
The effect on the selection efficiency, determined from events collected at random beam 
crossings, is a relative decrease of 5%. 

Finally, uncertainties in the simulation of the background were assessed by comparing 
the effect of each cut separately on the data and on the simulated backgrounds after the 
preselection. The relative differences, although compatible with the uncertainty due to 
limited statistics, were conservatively added in quadrature, for a total possible deviation 
of 9%. 

To derive the final result, all the above uncertainties were taken into account following 
the method of Ref . [35] . 

7.1.4 Results and interpretation 

The numbers of candidate events observed in the data between 1997 and 2000 are displayed 
in Table 3. These numbers are in agreement with the numbers of events expected from 
standard model background sources. 

In the framework of the MSSM with R-parity conservation, the outcome of this search 
can be translated into constraints in the (m^mg) plane when the stop quark decays with a 
100% branching fraction into a stable hadronizing gluino and a c quark. Regions excluded 
by this search at 95%C.L. are shown in Fig. 12, for 9 mix = 56° and 0°, corresponding to 
vanishing and maximal stop coupling to the Z, respectively. 
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Table 3: The numbers of candidate events observed in the data between 1997 and 2000, and the numbers 
of events expected from standard model background sources, for the small, intermediate and large AM 
selections. 



Year 


Luminosity [pb x ] 


Small AM 


Int. AM 


Large AM 


Obs. 


Exp. 


Obs. 


Exp. 


Obs 


Exp. 


2000 


207.3 


6 


8.1 


11 


14.7 


17 


17.0 


1999 


236.9 


9 


8.7 


15 


16.7 


19 


19.5 


1998 


173.6 


7 


6.9 


11 


13.4 


17 


15.2 


1997 


56.8 


1 


2.2 


5 


4.4 


4 


5.0 


Total 


674.6 


23 


25.9 


42 


49.2 


57 


56.7 



ALEPH eg 




stop mass (GeV/c ) 



Figure 12: The 95% C.L. excluded regions in the plane (mj, mg) from the acoplanar jet search at LEP 2, 
for maximal (horizontal hatching) and minimal (cross hatching) stop coupling to the Z. The shaded area 
corresponds to a stable stop and is not accessible to this search. 
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7.2 The AM ~ m D case 



The analyses of Sections 6 and 7. 1 do not suffice to efficiently cover the AM = m D line: 
along this line, the phase space for the t — > eg decay is so small that the stop lifetime may 
take any intermediate value, for which none of the previous selections is fully effective. 

Stop pair production followed by decays into eg was simulated for AM = as 
explained in Section 3, with stop proper decay lengths varying from f mm to f m. The 
acoplanar jet search of Section 7.1, the search for heavy stable charged particles of 
Section 6 and the search for kinks and secondary vertices of Ref. [3] were applied in 
turn to these simulated samples. Altogether, they allow all stop masses between 14 and 
80 GeV/ c 2 to be excluded for this very small mass difference, irrespective of the stop 
lifetime, for a vanishing stop coupling to the Z. 

8 Combined results 

As can be seen from Fig. 12, no absolute stop mass limit can be extracted from the results 
of the acoplanar jet plus missing energy search alone. An absolute mass limit is obtained 
by combining all searches presented in this paper in the following way. 

• Small gluino masses do not give rise to large enough missing energy to be addressed 
by the t — > eg search. However, the searches for e + e~ — > qqgg at LEP 1 presented in 
Ref. [17] and in Section 4 allow all gluino masses below 26.9 GeV/c 2 to be excluded 
when the gluino is the LSP (mg < mi). 

• Stable stops, i.e., with AM < tod, are excluded up to 15.7 GeV/ c 2 with the search 
for e + e~ — > qqqq at LEP 1 presented in Section 5, and up to 95 GeV/c 2 by the search 
for heavy stable charged particles presented in Section 6. 

• Promptly decaying heavy stops, i.e., with AM > mo, are excluded up to masses of 
85 GeV/c 2 by the acoplanar jet search of Section 7.1. 

• For AM ~ mo, masses up to 80 GeV/ c 2 are excluded independently of the stop 
lifetime, as mentioned in Section 7.2. 

The result of the combination is displayed in Fig. 13. All stop masses below 80 GeV/c 2 
are excluded at 95% C.L. when either the stop or the gluino is the lightest supersymmetric 
particle. When combined with the result of Ref. [2], all stop masses below 63 GeV/ c 2 are 
excluded irrespective of the nature of the LSP. 

Sbottom quarks were also searched for, but the acoplanar jet search proved not to be 
efficient enough to cope with (i) the four times smaller cross section of e + e~ ^bb — > bbgg, 
even when b-tagging is applied to reduce the background; and (ii) the very small mass 
differences (down to mx) for which the sbottom may decay promptly to sg. However, the 
Z lineshape limits and the excluded areas (1), (2) and (3) of Fig. 13 apply as well for all 
squark species, and in particular for the sbottom. All sbottom masses below 27.4 GeV/ c 2 
are therefore excluded at 95% C.L. when either the sbottom or the gluino is the lightest 
supersymmetric particle. 
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Figure 13: The 95%C.L. excluded regions in the plane (mj, mg) from the combination of all searches 
presented in this paper. The black area at very small gluino and squark masses is excluded by the precise 
measurement of the Z lineshape; Regions (1), (2) and (3) are excluded by the search for e + e~ — ► qqgg at 
LEP 1, for e + e~ — > qqqq at LEP 1, and for heavy stable charged particles from qq production at LEP 2, 
respectively. The hatched areas are excluded by the acoplanar jet with missing energy search at LEP 2. 



9 Conclusions 

Searches for stable, hadronizing squarks and gluinos have been performed with the data 
collected at centre-of-mass energies from 88 to 209 GeV by the ALEPH detector at LEP. 
No evidence for such a signal was observed in any of the production processes studied: 
e + e~ — > qqgg, e + e~ — > qqqq and e + e~ — ► qq. The following absolute mass limits were 
obtained at the 95% confidence level in the framework of the MSSM with R-parity 
conservation: 
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• a gluino LSP is excluded for mg < 26.9 GeV/c 2 ; 

• a down-type squark LSP is excluded for m, < 92 GeV/c 2 ; 

• an up-type squark LSP is excluded for m, < 95 GeV/c 2 ; 

• a sbottom quark NSLP is excluded up to masses of 27.4GeV/c 2 if the LSP is a 
gluino; 

• a stop quark NSLP is excluded up to masses of 80GeV/c 2 if the LSP is a gluino, 
and up to masses of 63 GeV/ c 2 irrespective of the nature of the LSP. 

The above squark and gluino LSP limits also apply in any supersymmetric model in 
which squarks or gluinos are long-lived. In particular, the scenario of Ref. [7] with a gluino 
of 12 to 16 GeV/c 2 decaying into a b quark and a long-lived sbottom with a mass of 2 to 
5 GeV/ c 2 is no longer viable. 

The results presented in this paper improve on related existing results [36, 39]. In 
particular, absolute lower limits on the masses of stable squarks, stable gluinos, and stop 
quarks decaying into stable gluinos, have been reported for the first time. 
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